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METHOD FOR IMAGE TEXTURE ANALYSIS 

FIELD OF THE INVENTION 

The invention relates to a computer-implemented method of analysing 
5 image texture. In particular, it provides a method of applying selected 

mappings to an image to generate a surface that may be analysed to 
characterise the texture of the image. 

BACKGROUND TO THE INVENTION 
10 image texture is a useful characteristic for object differentiation. 

However, image texture is a difficult characteristic to describe without using 
some kind of comparison to already familiar objects. Thus, human analysis 
of image texture is highly subjective. A computer based texture analysis 
method has the potential to provide a quantitative, and therefore more 
15 useful, measure of texture. 

Once a textured object is in digital form the problem of analysis 
becomes a mathematical problem rather than the verbal problem of 
description. The mathematical problem has proven to be no less. intractable. 
The problem of image analysis has traditionally been approached as 
20 an analysis of visual images or pictures of macroscopic objects, such as 

crops or forests, or microscopic objects such as stained tissues or cells. 
However, the principles of texture analysis are applicable to any image be 
it visual, acoustic or otherwise, image texture can be qualitatively evaluated 
as having one or more of the properties of fineness, coarseness, 
25 smoothness, granulation, randomness, Uneation or being mottled, irregular 

or hummocky. These are relative terms, a useful analysis of the 
characteristics of an image must include an absolute measurement of image 
properties. 

in general terms, the approaches to texture analysis can be divided 
30 into three primary techniques, these are : statistical, structural and spectral. 
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Statistical techniques include first order approaches such as autocorrelation 
functions and higher order approaches such as gray-tone spatial- 
dependence matrices (see Haraiick, Shanmugan & Dinstein; "Texturaf 
features for image classification", IEEE Transactions on Systems, Man and 
Cybernetics SMC-3(6), 610-621; 1973). Structural techniques include the 
use of primitives and spatial relationships (see Haraiick; "Statistical and 
structural approaches to texture", Proceedings of the IEEE 67(5), 78S-804; 
1 979). The main spectral technique is Fourier analysis. 

There are five steps in image processing. These steps are image 
acquisition, pre-processing, segmentation, representation and description, 
and recognition and interpretation. The invention falls into the step of 
representation and description. Persons skilled in the art will be aware of 
existing techniques for performing the remaining four steps of image 
processing. 

A general explanation of the prior art relevant to texture analysis can 
be found in Chapter 9 of "Vision in Man and Machine" by Levine, M; 
McGraw-Hill, Kew York (1985). Despite research in the area of texture 
analysis no single method has been developed that completely describes 
texture either verbally or mathematically. 

Abdel-Mottaleb makes mention of texture analysis in United States 
Patent No 5768333, as a technique for second stage discrimination of 
objects. While, Abdei-Mottaleb makes use of a large number of gray level 
thresholds, image texture is then characterised by computing many features 
of each connected component in the resulting threshold separately, no 2- 
dimensional surface is formed. 

OBJECT OF THE INVENTION 

It is an object of the present invention to provide a method of image 
texture analysis to usefully describe observed image texture. 

Further objects will be evident from the following description. 
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DISCLOSURE OF THE INVENTION 

In one form, although it need not be the only or indeed the broadest 
form, the invention resides in a method of analysing the texture of an image, 
including the steps of: 

mapping the image into a first set of binary representations of the 
image, wherein each representation in the set corresponds to a 
transformation by a first monotonically varying operator; 

mapping each binary representation of the first set of binary 
representations into a further set of binary image representations using a 
second monotonically varying operator, thereby forming a matrix of binary 
image representations; 

transforming the matrix by allocating a scalar value to each binary 
image representation; 

arranging the scalar values in a two dimensional space against 
parameters of the first monotonically varying operator and the second 
monotonically varying operator to form a first array of scalar values; and 

identifying defined texture characteristics from the resulting first array 
of scalar values. 

By monotonic operator is meant an operator acting on an image, the 
result of which only decreases (never increases) , with an increase in its 
parameter, or with a decrease in its argument. 

In preference, the first monotonicatiy varying operator is a threshold 
operator wherein the threshold operator applies a threshold value that varies 
monotonicaUy from a minimum value to a maximum value to produce the f irst 
set of binary representations. 

In preference, the second monotonically varying operator is a spatial 
operator wherein the spatial operator applies an opening that varies 
monotonically from a minimum value to a maximum value to produce the 
further set of binary representations.' Suitably, the second monotonically 
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varying operator is a morphological opening that removes ail image objects 
below a defined size. 

The method may include the further step of calculating a second array 
of scalar values derived from the first array of scalar values by taking at each 
point the minimum of the first difference in a first direction and the first 
difference in a second direction. 

The step of identifying defined texture characteristics may include the 
further step of comparing the first or second array of scalar values to a 
knowledge base of first or second array values known to indicate defined 
texture characteristics. 

The method may include the further step of calculating one or several 
texture feature values by computing linear sums of the values in the first or 
second arrays of scalar values. In preference, the coefficients of these linear 
sums are determined by a principal components analysis on the data. 

In a furiherform, the invention resides in an image analysis apparatus 
comprising: 

storage means for storing a digital representation of an image; 

processing means for : 

mapping the image into a first set of binary representations of 
the image, wherein each representation in the set corresponds 
to a transformation by a first monotonicaliy varying operator, 
mapping each binary representation of the first set of binary 
representations into a further set of binary image 
representations using a second monotonicaliy varying 
operator, thereby forming a matrix of binary image 
representations; 

transforming the matrix by allocating a scalar value to each 
binary image representation; 

plotting the scalar values in a two dimensional space against 
parameters of the first monotonicaliy varying operator and the 
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second monotonically varying operator to form a first array of 
scalar values; and 

identifying defined texture characteristics from the resulting 

first array of scalar values; and 
display means for displaying the first array of scalar values and 
identified texture characteristics. 

BRIEF DETAILS OF THE DRAWINGS 

To assist in understanding the invention preferred embodiments will 
now be described with reference to the following figures in which : 
FIG 1 is a flow chart showing the method of image texture analysis; 

FIG 2 shows the array of values of the granotd surface; 

FIG 3 displays a granold surface; 

FIG 4 is an example of the application of the method to an image; 

FIG 5a shows the granold surface derived from the example of FIG 4; 

FIG 5b shows the granold spectrum derived from the granold surface 

of FIG 5a; 

FIG 6 shows a first example of the application of the method to an 

image of ammonium sulphate pellets; 
FIG 7 shows a second example of the application of the method to an 

image of ammonium sulphate pellets; 
FIG 8 shows a third example of the application of the method to an 

image of ammonium sulphate pellets; and 
FIG 9 is a block diagram of an apparatus for working the invention. 

DETAILED DESCRIPTION OF THE DRAWINGS 

In the drawings, like reference numerals refer to like parts. Referring 
to FIG 1, the flow diagram indicates the steps of a method of texture 
analysis. The method is preferably implemented in a computer, persons 
skilled in programming will be able to readily produce appropriate code in a 
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suitable language to implement the method. 

The flow diagram is directed to the situation in which the monotonic 
mappings are gray ievel thresholding and spatial opening. The example is 
two-dimensional but the method is applicable to multi-dimensional 
applications. A mathematical description is presented later which generalises 
the specific example to a wider range of situations. 

In the fiow diagram, the gray level image to be analysed is designated 
as X. A threshold operator <2> i5 is applied at a first value of I = i min . Typically, 
an 8-bit gray level image might first be thresholded at gray level 0. A binary 
image 0,(X) is obtained by this thresholding operation where pixels having 
a gray ievel above the threshold will have a value of one and all other pixels 
have a value of zero. Applying spatial operator x ¥ j then morphologically 
opens the binary image. 

In the example described here, a set of disc shaped structuring 
elements of increasing radius is used as the morphological opening. Thus, 
the smallest value of j = j mjn corresponds to the first member In the set, which 
has the smallest radius. All pixel clusters having a radius less than j are set 
to zero and the rest are retained. The remaining area of the resulting binary 
image ^-^(X) is then calculated (for example, by counting all pixels having 
a value of one) to give a value e^O^X). This value is stored in a two- 
dimensional array G x in a position corresponding to the I th threshold operator 
versus the spatial operator. 

The value of j is iterated and the next spatial operator x F j+1 is applied 
to <£>i(X) and the value S^^O^X) is calculated and stored in the array at 
location G x (i,j+1). The iteration and application of the spatial operator 
continues until j = j m3X . The value of j max may be arbitrarily chosen or may be 
the value of j such that dW^^X) - 0. This process of continually opening 
a binary image with a suitable set of structuring elements and measuring the 
remaining area produces a morphological granulometry of the binary image. 

Once the granulometry of the first binary image has been completed 
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the original gray level image is thresholded for the next value of 1. The 
morphological granulometry for the binary image is then formed and 

entered into the granoid surface array. The process is repeated, as shown 
in FIG 1 , until ! = i max . The value for I may be chosen arbitrarily or may be the 
value of the highest gray level threshold possible in the image. For an 8-bit 
gray level image this would be 255. 

The resultant array is represented in FIG 2 and is referred to as the 
topographical granoid surface. The array has (i ma5f - i m(n ) x G max -j mSn ) elements 
that can be plotted as a surface in the manner shown in FIG 3. Corners 
appear on the surface that indicate a relationship between size and gray 
level of granules in the original image. An idealised corner is shown in the 
surface G x (i,j) of FIG 3 at a position corresponding to l=i c and j=j c . Further 
processing, for example by taking at each point the minimum of the first 
difference in the I direction and the first difference in the j direction of the 
granoid surface, can be performed to form the granoid spectrum (also called 
the granoid korner surface). As shown in the following figures, further 
processing can result in a spectral-like result that simultaneously reveals the 
distribution of size and gray level of granules within the original image. 

FIG 4 shows an example of a number of images produced by the 
method described above. FIG 4a shows a 50x50 pixel gray level image. FIG 
4b is a binary image produced by thresholding the image of FIG 4a at a 
threshold level of 172. The image of FIG 4b is opened with a disc shaped 
structuring element of radius 1 to produce the image of FIG 4c. The image 
of FIG 4c is opened with a disc shaped structuring element of radius 2 to 
produce the image of FIG 4d. This process is continued with increasing 
radius until the remaining image is reduced to zero. The original image is 
thresholded at a gray level of 209 to produce the binary image of FIG 4e. 
Disc shaped structuring elements of radius 1 and 2 are used to produce the 
images of FIG 4f and FIG 4g respectively. 

FIG 4 clearly demonstrates that the number of remaining pixels 
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decreases with increasing threshold and increasing radius of the spatial 
operator. The actual values of the granold surface calculated for each image 
are G x (1 72, 1 )=1296, G x (1 72,2)=1210, G x (209,1 )=1 33 and G x (209,2}=0. 

The complete granoid surface G x (i,j) for the image of FIG 4a is 
displayed in FIG 5a. The surface was constructed using a threshold operator 
that applied gray ievei thresholds from 0 to 255 and a spatial operator that 
applied a set of disc shaped structuring elements with radii from 1 to 1 6. The 
plot has been normalised to the original image area of 2500 pixels. The 
corresponding granold spectrum K^ij) is shown in FIG 5b. The granold 
spectrum provides a spectra!-!ike plot that characterises the original image 
in a reproducible and objective manner. 

EXAMPLE 

Characterisation of matter by pellet size and colour intensity may be 
applicable to a variety of industries, such as confectionery, pharmaceuticals 
and fertilisers, for quality control. FIG's 6-8 provide an example of the 
application of the invention to the characterisation of ammonium sulphate 
pellets. 

In each of FIG's 6-8, image (a) is a gray level image produced by a 
video camera and frame grabber as ammonium sulphate pellets fall from a 
drum granuiator in a fertiliser plant The pellets of F1G6 and FIG 7 are the 
same size but the pellets of FIG 8 are larger. Green food dye was added to 
the pellets to provide colour variation with FIG 6 being the darkest and FIG 
8 being the lightest. 

Image (b) of each figure is a plot of the granold surface produced by 
applying a threshold operator that applied gray level thresholds from 0 to 255 
and a spatial operator that applied a set of 16 disc shaped structuring 
elements. Image (c) is the corresponding granold spectrum. 

Further processing of the granold spectrum involved identifying the 
threshold/radius location of the maximum value of the surface. The value 



WO 00/52647 



PCT/AUOO/00141 



9 

indicates the gray level and size of the majority of the pellets in the image. 
!t wiii be appreciated that statistical analysis of the granoid spectrum can 
characterise the pellets by the mean and standard deviation of the size and 
colour of the pellets. Such an analysis provides a useful standard for quality 
control. 

It will be appreciated that the foregoing examples are limited to two- 
dimensional space but can be generalised to multi-dimensional images to 
produce a similar granoid surface and granoid spectrum. Furthermore, the 
operators that produce the mappings are not limited to threshold, opening 
size and area. Other suitable operators may include contrast dynamic, 
attribute opening, etc. 

The following generalised description is presented mathematically in 
a manner that will be understood by those skilled in the relevant art. 

MATHEMATICAL DESCRIPTION 

In order to describe the method it is necessary to define certain 
notations. In the following, "R + " denotes the non-negative real numbers 
(including zero), "Z" the integers {...,-1,0,1,...} and "N" the natural numbers 
{0,1,2,...}. 

We define 

{F}:2L" -> Z» (1) 
Here {F} represents the general class of m-integer valued functions on n- 
dimensionai space. We would typically use such a function to represent a 
gray scale 2-D image, in which case we would write, for example, "let 
X e F be a 2-D gray scale image This statement implies that n=2 and 
/77=1 . We keep the dimensions of the image representation general (at n and 
m) so as not to unnecessarily limit the scope of the method to 2-D images. 
The domain of the image X e F will be denoted by Dom(X) . The domain 

is the set of al points in Z"for which X is defined. 
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Next we define, 

B e />(Z") (2 ) 
Where P is the powerset operator, so ?(Z" ) is the set of al! subsets of Z" . 

B represents the general class of binary images on ^-dimensional integer 
space. We would typically use such a notation to represent a 2-D binary 
image where we would write, "let be a 2-D binary image ..." which statement 
implies that n+2. Again, we keep the dimensions of the image general at n 
so as not to unnecessarily iimit the scope of the method to 2-D images. The 
domain of the binary image 5eB will be denoted JDomQ i?) . 

Based on the image spaces defined above we now define the 
mappings or operators used to form the granold surface. Firstly, we define, 

O :JF -» B (3) 
This function to set mapping takes a gray scale image and returns a binary 
image. As an example let JeFbea 2-D gray scale image, so we could 
write, B = <D {X ) . The resulting image B e B is a 2-D binary image. 

The thresholding of regions based on region minima or maxima 
values, the thresholding of regions based on contrast dynamic values, and 
the gray level threshold operation, well known in image processing, are all 
examples of this type of mapping. 

Secondly we define, 

^:B -» B (4) 
This set to set mapping, which we will call an "operator", takes a binary 
image and returns a binary image. As an example, let B e B be a 2-D binary 
image so we could write, A = *P (B) . The resulting image A e B is also a 
2-D binary image. 

The binary AND operation, the attribute opening operation, and the 
morphological binary opening operation, also well known in image 



WO 00/52647 



PCT/AUO0/O0I41 



11 

processing, are all examples of this type of mapping. 

It is the intention that these mappings preserve the domain or support 
of the images in ail cases, that is, 

Dom(& (X)) = DomQX) (5) 

5 and 

DomQ¥ (B)) = Dom{B) (6) 
and with the composition of these mappings' 

DomQ¥ O (X)) = Dom(JsT) (7) 
Lastly, we define the set measure, 
10 £:B-> R + (8) 

This set to scalar mapping, which we will cail a "measure", takes a binary 
image and returns a non-negative scalar value. The area (number of object 
pixels) of a binary image is an example of this type of mapping, so we could 
write, for example, area = $(B) . 

15 The next step is to parameterise the mappings <& and ^ to give 

families of mappings and operators Indexed on non-negative integer 
parameters i e N . Therefore we have families, 

{O,}:^-* B (9) 

and, 

20 {^}:B -» B (10) 

It is important for the method that the above mappings posses 
additional mathematical properties as follows: 

Property 1 - Monotone Decreasing for Parameter. 
25 The families of mappings O, and X F J must be "monotone decreasing" 

with respect to their parameter I. Also, to ensure that (in the limit) the 
mappings extract all the information from the images, there should exist 
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values M and N (which may depend on the images) such that for i>M and 
i>N respectively, the appropriate mappings produce the empty-set, 0 . In 
contrast, we should afso ensure that the index zero mappings have the least 
possible effect on the image and respectively produce the domain of the 
image or reproduces the image itself. 

in mathematical terms, for the mappings d> f and for any X e F , 



i = 0 => $ . (X) = Dom(X) (l l) 

0<i<j<M => (D/AOc $,(^) (12) 

i> M => O i (^)= 0 (13) 

Likewise, for the operators and for any ^eB, 

/=0 =>W i tA)=A (14) 

0<i<j<N =>^(^)s^,(^) (15) 

/> AT ¥,(-4) = ^ (16) 



Property 2 - Monotone Increasing for Argument 

The families of operators (set to set mappings) T,. need to have some 
additional properties to ensure the granold surface is well behaved. Firstly, 
V F J needs to be "anti-extensive" and to remove detail from the binary image 
B e B , that is, for all i € N 

^XB)=B (17) 

Secondly, the operators ^ . must be "monotone increasing" with 
respect to their arguments, the binary images, that is for the binary images 
A,B & B for a fixed integer i e N 

A<= B=> ¥ ; .(^)c VXB) (18) 
Property 3 - Measurement 
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Finally, the measure S , to represent the usual sense of measuring 
something, must be : (I) zero if the object to be measured is empty, 
£(0)=O (19) 
and (ii) increasing, that is, 

A c= B^> S(A)< S(B) ■ (20) 

Now we consider the composition of these mappings as follows 5 
(^O,):^-* R + (21) 

This function to scalar mapping, which we have previously termed a 
"granold" (granulometry of threshold), takes a gray scale image and returns 
a non-negative scaiar which depends on the values of the parameters 7 and 
j (and of course the image itself). 

To show the dependence on I and j explicitly, we choose to denote 
the granoid (z',y):N 2 -» R + of an image X € F as a 2-D function, 

def 

G x Q,j) = {cW J <t>,XX) (22) 

If G is represented as a function on two dimensions or a surface then 
regions I and j both have a large effect on the value of 6 wili appear as 
"corners" in the surface. The location and "sharpness" of these corners will 
reveal the particular values of / and y that have a large effect on the value of 
6. 

The following first differences are useful for analysing G. These 
differences approximate the directional derivatives of the granold and will be 
large wherever there is a rapid fall in the granold in the corresponding 
directions. 

def 

h' x QJ) = G x (i,fi- G x (i+ 1,7) (23) 

def 

& y x (i,j) = L> x Q,n-G x Qj+Xy (24) 
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The quantity, 

def 

K x O,j) = min(A x A ,(/,y),A^(/,y)) (25) 

is useful to determine where the granoid surface has simultaneous large first 
differences in both the d> ; and the T ; directions, and produces the granoid 
spectrum. 

To classify textures we may compute the texture features, 

S k = £ a fa/ G>(z, j) (26) 

if 

and/or, 

T* = HPuffKxQ.J) (27) 

if 

15 where the coefficients, a k8 and B kjj and would be found by, for example, 

principal components analysis of the data. 

A suitable apparatus for putting the invention into effect is shown in 
Fig 9. An image is acquired with acquisition means 1 . This could be a digital 
camera or digitizer that converts an analogue signal. Pre-processing may 

20 occur at 2. This may be on-board processing in the acquisition means 1 , or 

may be, for example, post-acquisition filtering. 

Signals are passed to a processor 3 that is programmed to work the 
method described above. Required parameters may be stored in ROM 4 and 
temporary storage is provided by RAM 5. The result of the image texture 

25 analysis is displayed on monitor 6. Other data storage means such as disk 

or tape (not shown) may also be connected to the processor 3. The data 
storage may store a knowledge database of previous image texture analyses 
that can be used to identify texture characteristics by comparison. 

Persons skilled in the relevant art will appreciate that the method of 
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image texture analysis can be applied to a diverse range of situations. The 
inventors have demonstrated the application of the method to fertiliser 
pellets. Preliminary considerations indicate that it is also applicable to 
analysis of cervical celis on a Pap smear slide on a microscopic level and 
5 skin screening at a macroscopic level. Throughout the specification the aim 

has been to describe the preferred embodiments of the invention without 
limiting the invention to any one embodiment or specific collection of 
features. 
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